Structure of the Local-field factor of the 2-D electron fiuid. Possible evidence for 

correlated scattering of electron pairs. 
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(Dated: 2-june-2003) 

The static local-field factor (LFF) of the 2-D electron fluid is calculated nonperturbatively using 
a mapping to a classical Coulomb fluid [Phys. Rev. Lett., 87, 206404 (2001)]. The LFF for the 
paramagnetic fluid differs markedly from perturbation theory where a peak near 2kF is expected. 
Our LFF has a quasi-linear small-k region leading to a peak close to Skp , in agreement with available 
quantum Monte Carlo data. The structure in the LFF and its dependence on the density and 
temperature are interpretted as a signature of correlated scattering of electron pairs of opposite 
spin. The lack of structure at 2A;f implies weakened Friedel oscillations, Kohn anomalies etc. 

PACS numbers; PACS Numbers: 05.30.Fk, 71.10.-|-x, 71.45. Gm 



Introduction. — The physics of the uniform two- 
dimensional electron fluid (2DEF) depends crucially on 
the "coupling parameter" F — (potential energy) /(kinetic 
energy). The F for the 2DEF at T = and mean den- 
sity n is equal to the mean-disk radius rg — (Trn)"^/^ per 
electron. The parameter Vg, the spin polarization ^ and 
the temperature T are the only variables in this problem. 

The response function xi^^^^) is a property of the 
2DEF sensitive to exchange-correlation effects. It is ex- 
pressed in terms of a reference "zeroth-order" 
and a local-field factor (LFF), denoted by G{k,oj). 

X(fc, u;) = x%{k, u:)/[l - Vk{l - G{k, u;)}xl{k, u;)] (1) 

The LFF is closely related to the vertex function A(fc, oj). 
The static form, G{k), is identical with G{k,0). As 
such, considerable effort has been devoted to determining 
G'(fc), using perturbation theory, kinetic-equation meth- 
ods 0,01, etc. A partially analytic, semi-empirical ap- 
proach invokes parametrized models constrained to sat- 
isfy sum rules isj which fit (0, IE 0) to Hmited re- 
sults obtained from quantum Monte-Carlo (QMC) sim- 
ulations While these efforts have considerably ex- 
tended the available data, it is still restricted to the fit- 
ted Ts regime. Also, these fits usually do not invoke any 
physical model. This truly emphasizes the difficulty and 
delicateness involved in the determination of the LFF. 

In this paper we show that the classical LFFs, eval- 
uated for a classical Coulomb fluid which is an approx- 
imate mapping of the quantum fluid, agree remarkably 
well with the available quantum data, reproducing the 
quasi-linear behaviour in the k region up to about 2 kp, 
which gets further extended into a peak structure near 
3kp. 

The classical mapp ing was discussed in a number of 
papers 0, 0, 0, 0| where we showed that the static 
properties of the 2D and 3D electron systems, (or even 
electron-proton systems 0|), can be calculated via an 
equivalent classical Coulomb fluid having a temperature 
Tq such that it has the same correlation energy Ec as the 



quantum system at the physical temperature T = 0. The 
mapping is based on an extension of the classical Kohn- 
Sham equation at Tq so as to mimic the quantum sys- 
tem. The "quantum temperature" Tq in 2-D was found 
to be [13, 



t = Tq/Ep = 2/[l + 0.86413(: 



,1/6 



1)^ 



(2) 



where Ep — l/r^ is the Fermi energy in Hartrees. At 
finite T, the classical-fluid temperature T^f is taken 
to be (T2 -f- T2)1/2^ as discussed in Ref. 0]. The 
pair-distribution functions (PDFs) are given by the 
hyper-netted-chain (HNC) equation inclusive of 

bridge terms. The HNC equations and the Ornstein- 
Zernike(OZ) relations are \\^ '. 

gij{r) = exp[-Pcf(t)ij{r) + h,j{r) - c,j{r) + Bij{r)] 

hij{r) = c.y(r) + SjTis / dr'hi^s{\r - v'\)cs^j{r') (3) 



These involve: (i) the pair-potential n^ij \ 



j{r), (ii) the bridge 

function -Bjj (j') [ij, [ISj- The other term, e.g, c{r), is 
the "direct correlation function" . These are discussed in 
ref. |lOllTl| . and briefly below. This classical mapping of 
quantum fluids within the HNC was named the CHNC. 

In effect, although S{k) and related properties (e.g, 
g{r), LFF) of a quantum fluid have to be determined 
(traditionally) by first evaluating S{kui), and then in- 
tegrating over w to obtain S{k), the CHNC mapping 
leads directly to good estimates of S{k), g{r) etc. In 
this paper we show numerically that the classical LFFs 
obtained from CHNC are in remarkable agreement with 
QMC data for the available values. This confirms our 
basic premise that the static properties of the classical 
fiuid provide a good approximation to the corresponding 
properties of the quantum fluid. 

The 2-D LFFs do not have the form indicated by 
standard perturbation theory 6J. Such calculations give 
LFFs with a "hump" at 2kp . Here we find that the inter- 
actions have extended the quasi-linear region and pushed 
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the usual 2kp hump towards ^ 3kp. The hmited set of 
QMC data for the 2-D LFF Q agree with this. We exam- 
ine the behaviour of the hump in the LFF, with and with 
out the clustering term (bridge term), and as a function 
of and T, and find that the hump in the classical LFF 
results from mp-down electron correlations. This sug- 
gests that correlated scattering of singlet pairs may play 
a role in the quantum fluid as well, since the usual 2kp 
anomaly arises from scattering of uncorrelated electrons 
across the Fermi-disk. 

The local-field factor. — We consider the static form 
G(fc), defined with respect to a reference "zeroth-order" 
response function. The Lindhard function Xl(^) ^ 
ten used for this purpose. However, another natural 
choice 0, 0| is to use the "density-functional" non- 
interacting form xC^)" containing the occupation num- 
bers corresponding to the interacting density. The two 
choices mainly affect the large-fc behaviour of the LFF 
01 . Thus, for the Lindhard reference used in QMC, 



lim G(fc) 

k — >oo 

e'en 



Coofc + l-.g(0) (4) 
-ir,a/2)d[r,e,ir,)]/dr,. (5) 



We showed in Ref. that the LFFs from CHNC have 
the 1 - g{0) hmit. The CHNC data and the QMC 
data can be compared on the same footing by remov- 
ing the asymptotic linear-fc dependence. Here edTs) is in 
Hartrees/electron, k is in units of kp, and a — l/^/{2). 
The remaining term, 1 — g(0) depends on the estimated 
"on-top" value of the PDF, a subject discussed by Bu- 
lutay et al. ^3 the context of the CHNC and other 
methods. 

The CHNC provides a very simple formula for the LFF, 
via the classical-fluid. Unlike in the quantum case, for a 
classical fluid, is directly related to the structure 

factor. 



S,,(fc)=-(l//3)x,,(fc)/(n,n,)i/2. 
Hence, taking the paramagnetic case for simplicity, 



Vc{k)G{k) = Vc{k) - ^ 
n 



1 



1 



S{k) S"{k) 



(6) 



(7) 



Here Tcf equals Tq if the physical temperature T = 0. 
In CHNC the x°(fc) and S'°(fc) are based on a Slater 
determinant and not on the Lindhard function. QMC 
results use a reference Xl such that the LFF contain a 
kinetic-energy tail, as discussed in Eq. 0] 

The S{k) needed in Eq. 0is explicitly known from the 
CHNC calculation. Alternatively, any other source of 
e.g., QMC, may be used, while S°ik) for the 2DEG 
is analytically known. The Coulomb potential Vc{r) oc- 
curring in Eq. \7\ needs explanation. The Coulomb op- 
erator for point-charge electrons is 1/r. However, the 
classical electron at the temperature T^f is localized to 
within a thermal de Broglie wavelength. Hence, the ef- 
fective classical interaction in CHNC is the "diffraction 



corrected" form 

V^r) = (l/r)[l 



y,(fc) = 2^[fc-i-(fc,\ + fc^)-V2] 



(8) 
(9) 



By numerically solving the Schrodinger equation for a 
pair of 2-D electrons in the potential 1 /r and calculating 
the electron density in each normalized state j '^^ S^t 

where Tcf is in au. Here fc^^ is the familiar 3-D form 
of the de Broglie wavevector, (27rm*rg)^/^, where m* is 
the effective mass of the electron pair. We emphasize that 
the G{k) can be calculated only if this modified potential 
Vc{k) were used. The large-fc sum rule for G{k) is recov- 
ered with this Vc{k) which also contains the Tcf entering 
into Eq. 13 In effect, the classical equation successfully 
satisfies the small and large-fc sum rules satisfied by the 
quantum G(fc) itself. 

The Coulomb potential Vc becomes large for small fc, 
and explicit cancellation of Vc by the terms in the 1/5 — 
l/5o is desirable in numerical work. This can be done 
by rewriting the structure factors in terms of the direct 
correlation functions using the OZ relations [T^. 

The bridge term models short-ranged correlations (fc > 
2kp) in the 2-D LFF. Hard-sphere (0) or hard-disk 
models can be used to obtain an explicit form for the 
bridge function 10, 15]. This gives satisfactory results, 
esp. for strong coupling where it is most needed. How- 
ever, unlike in 3-D classical fluid-studies have 
not used the bridge function to fit the compressibility 
sum rule. The latter is approximately satisfied even with- 
out the bridge term. It was used as a short ranged inter- 
action to mimic the back-flow effects [13 of QMC, and 
comes into play mainly for fc > 2kp. An improved bridge 
function would be based directly on the Vc(r) rather than 
the hard-disk model used here. The required hard disc 
diameter a = 2rsy^r], where rj is the packing fraction. It 
is given by 



?7 = 0.235ry V[l + 0.86413(r]/^ - 1) 



(10) 



and is based on the dependence required by the Gibbs- 
Bogoliubov form of the free energy with respect to a ref- 
erence fluid. The bridge term becomes zero when 77 is set 
to zero, and hence we can study the LFF with and with- 
out cluster terms. Also, since the Vc{r) becomes small 
for r ^ 0, the large-fc behaviour of the LFF based on the 
r] — bridge function would be of interest. Similarly, for 
small-fc (i.e, large-r), the short-ranged correlations are 
irrelevant and the = behaviour is retained. 

Comparison of the CHNC results with QMC data. — 
In Fig. 1 we compare the CHNC results for rg = 2, 5, 
and 10 (with and without the bridge term), with the 
QMC results. The CHNC results for the LFF have not 
been fitted to any outside data. In comparing with QMC, 
it is necessary to remove the large-fc dependence arising 
from the Lindhard reference (see Eq. 0)1 . The subtrac- 
tion of the Coofc term must be applied asymptotically 
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and this leads to some arbitrariness in deciding on the 
"asymptotic" regime. Atwal et al. |^ use what they call 
"an admittedly ad hoc" scheme in their Eq. (22). In 
Fig. 1 we give the QMC data (labeled 'qmc-a') where the 
asymptotic term of Davoudi et al. Q has been used as 
a means of removing the /c-asymptote. It modifies even 
the data points smaller than 2kp. After subtraction, the 
maximum-amplitude LFF occurs for Ts = 10, rather than 
at = 5.41, as was the case prior to subtraction. 

The LFF from CHNC without the bridge term {tj = 0) 
shows a large quasi-linear behaviour of the LFF for k up 
to and beyond 2k p. The bridge term extends the quasi- 
linear region and introduces short-range effects (i.e, for 
larger-A:) and producing a hump, agreeing with QMC 
data, even though the QMC /c-range is limited. The 
CHNC data with and without bridge terms go to the 
CHNC 1 -g(0) limit for large k. The CHNC-bridge LFF 
has oscillations in the pre-asymptotic beyond-the-peak 
region. This is not seen in the QMC points. QMC seems 
to follow the 77 = curve for large-/c. This is probably 
realistic since Vc{r) near r — s- is much softer than the 
hard-disk potential used by us for the bridge term. 

Thus we see that the classically calculated LFF of the 
CHNC Coulomb fluid provides a good approximation to 
the QMC generated LFF in all the available cases. Hence 
we can use the CHNC to calculate classical LFFs for 
other Tg values where QMC data are not available. 

The featureless 2kF region and the appearance of a 
broad hump near Skp are shown for rs=5,10,20, and 30 
in the lower right panel. The increased coupling (larger 
Ts) moves the peak to s/iort er wavevectors. The absence 
of structure near 2k p implies a weakening of Friedel os- 
cillations and Kohn anomalies in 2D. 

Discussion — The strong coupling effects in the 2DEG 
were modeled in CHNC via a bridge term limited to 
short-range interactions ( k > 2k p). This term plays a 
role only for anti-parallel spins, i.e., in gi2{r). Also, avail- 
able results (not discussed here) show that the hump near 
Skp does not appear for the parallel spin case (the anti- 
symmetric LFF), where the 2kF behaviour is similar to 
that expected from low-order perturbation theory. These 
considerations suggest that the structure near Skp my be 
a result of correlated pair processes. The broadness of the 
peak suggests that this is not a sharp process. The lack 
of Pauli exclusion between two opposite-spin electrons 
and the strong coupling would lead to correlated pairs. 
Our results show their importance in the classical fluid 
which is the CHNC map of the quantum fluid. Hence 
we make the hypothesis that such pairs play a role in the 
quantum liquid as well. In Fig. |2Ia) an uncorrelated elec- 
tron scatters with an electron across the Fermi disk, with 
a change of momentum Afc = 2kp. The structure seen 
in the LFFs near 2kp in weak coupling arises from such 
uncorrelated scattering across the Fermi disk. Although 
there is another electron of opposite spin in the same kp- 
state, it is uncorrelated with the scattering electron and 
takes no part in the transition. 

Consider the correlated case, Fig. |21(b), where the up- 
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FIG. 1: The local-field factors calculated from CHNC, with 
and without (r; = 0) the bridge term. The QMC data, 
squares, ref Q were extracted from Davoudi et al.^^, with 
the asymptotic-fc behaviour subtracted out (see text), and la- 
beled "qmc-a". The lower right panel shows the ^ S/cf hump 
structure of the LFF for rs=5,10,20 and 30 obtained from 
CHNC. 




(b) (a) 



FIG. 2: Electron scattering across the Fermi disk, (a) Single- 
electron scattering, (b) Scattering of a correlated u- and d- 
spin pair. (c)A typical Feynman graph describing these cor- 
relations. 

spin u, and down-spin d electron are in two states at kp^ 
making an angle a in the Fermi disk. If the Coulomb re- 
pulsion were absent, the u, d pair would occupy the same 
state with a = 0. Unlike in the uncorrelated case (a), 
scattering of the correlated u, d pair can occur in a con- 
certed manner and would lead to a Afc — Akpcos{a/2). 
Depending on the correlations existing in the fluid, some 
optimal value of a would be most probable. The max- 
ima in the LFF for rs=5, and 20 are at Ak/kp= 3.34 
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FIG. 3: Temperature {t = T/Tp) dependence of the LFF 
peak at rs=20. 

and ~ 3. This mimics the shift in In the peak in the 
S{k) from 2.98 to 2.68 for rs= 5 and 20. If these pro- 
cesses are to be treated in a diagrammatic theory for the 
polarization operator, then diagrams hke Fig. |21 (c) are 
needed. The peak structure is strongest at the density 
Tg = 5.41, where correlated singlet-pair effects are most 
strong. After that exchange effects begin to counter cor- 
relation effects. 

At finite temperature, correlated-pair scattering 
should become weaker and the peak position should shift 



towards larger k/kp, rather than towards 2fci?. This is 
confirmed in the finite-T data for rs=20 shown in Fig. El 

Finally we remark that the bridge term ( similarly, the 
back-flow term in QMC) provides extra pair-interactions 
which make the paramagnetic fluid energetically less 
favourable than the ferromagnetic phase. Thus the re- 
cently proposed para— *■ ferro transition 0, |23| occurs 
only if bridge contributions (or, in the QMC case, back- 
flow terms) are included in the analysis. That is, there 
are no para-^ ferro transitions in the rj = CHNC cal- 
culation [r^ . 

Conclusion — We have shown that the CHNC derived 
LFF provides a remarkably good representation of the 
quantum simulation results so far available. Unusual fea- 
tures of the 2D-LFF not found in the 3D-case, and unex- 
pected from perturbation theory, were examined via the 
CHNC method. The lack of structure near 2kp and the 
presence of unexpected structure near Skp which arises 
only when cluster-terms are included in the classical map 
suggest them to be signatures of correlated singlet-pair 
scattering in the 2-D electron fluid. The possibility of 
such scattering would be very relevant to theories of su- 
perconductivity in 2-D systems, spintronics and related 
topics. The CHNC method thus provides a useful ex- 
ploratory tool for strongly correlated regimes inaccessi- 
ble by standard analytical methods. On-line access to 
our CHNC codes and more details may be obtained at 
our website[2lj. 



[t] electronic mail address: chandre@nrcphyl.phy.nrc.ca 

"*NRC Visiting scientist. 
[1] A.K Rajagopal and J. C. Kimball, Phys. Rev. B 15, 2819 

(1977) ' 

[2] K. S. Singwi,M. R Tosi, R. H. Land, and A. Sjolander, 

Phys. Rev. 176 589 (1968); 
[3] N. Iwamoto, Phys. Rev. A 30, 3289 (1984) 
[4] B. Davoudi, M. Polini, G. F. Giuliani, and M. P. Tosi 

Phys. Rev. B 64, 153101 (2001) 
[5] J. Moreno and D. C. Marinescu, corid-mat/0206465 
[6] G. S. Atwal, L G. Khalil and N. W. Ashcroft, Phys. Rev. 

B, 67 115107 (2003); I. G. Khalil, M. Teter and N. W. 

Ashcroft, Phys. Rev. B 65, 195309 (2002) 
[7] S. Moroni, D. M. Ceperley, and G. Senatore, as quoted 

in Ref. 10 of B, Davoudi et al., Ref. Q. 
[8] M. W. C. Dharma-wardana and F. Perrot, Phys. Rev. 

Lett. 84, 959 (2000) 
[9] Frangois Perrot and M. W. C. Dharma-wardana, Phys. 

Rev. B, 62, 16536 (2000) 
[10] Frangois Perrot and M. W. C. Dharma-wardana, Phys. 

Rev. Lett. 87, 206404 (2001) 
[11] M. W. C. Dharma-wardana and F. Perrot, Phys. Rev. 

Lett. 90, 136601 (2003) 
[12] M. W. C. Dharma-wardana and F. Perrot, Phys. Rev. B, 

66, 14110 (2002) 
[13] J. M. J. van Leeuwen, J. Groneveld, J. de Boer, Physica 



25, 792 (1959) 

[14] F. Lado, S. M. Foiles and N. W. Ashcroft, Phys. Rev. A 

26, 2374 (1983) 

[15] Y. Rosenfeld, Phys.Rev. A 42, 5978 (1990), M. Bans et 
al., J. Phys. C: Solid State Phys. 19, L463 (1986) 

[16] G. Nicklasson, Phys. Rev. B 10, 3052 (1974), A. Holas, in 
Strongly Coupled Plasma Physics, Edited by F. J. Rogers 
and H. E. DeWitt (Plenum, New York,1987) p.463, G. 
Vignale and K. S. Singwi, Phys. Rev. B 32, 2156 (1985) 

[17] C. Bulutay and B. Tanatar, Phys. Rev. B 65, 195116 
(2002) 

[18] M. W. C. Dharma-wardana and F. Perrot, cond- 
mat/0211127 This paper contains preliminary results 
from somewhat less accurate numerical procedures, 
partly associated with ensuring that the bridge function 
was truly short ranged. We find that the best way to en- 
sure a short-ranged bridge function is to use no bridge 
function {rj = 0) up to fc 2fe_F , and then include a finite 
rj for k > 2kF- The cross over is quite smooth, esp. at 
large coupling. 

[19] Y.Kwon et al., Phys. Rev. B 48, 12037 (1993) 
B. Tanatar et al., Phys. Rev B 39, 5005 (1989) 

[20] Attacalite et al.,Phys. Rev. Lett. 88, 256601 (2002) 

[21] http://nrcphyl.phy.nrc.ca/ims/qp/chandre/chnc, 



